Kondo Correlations and Fano Effect in Closed AB-Interferometers 
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We study the Fano-Kondo effect in a closed Aharonov-Bohm (AB) interferometer which contains 
a single-level quantum dot and predict a frequency doubling of the AB oscillations as a signature of 
Kondo-correlated states. Using Keldysh formalism, Friedel sum rule and Numerical Renormalization 
Group, we calculate the exact zero-temperature linear conductance G as a function of AB phase ip 
and level position e. In the unitary limit, G(</s) reaches its maximum 2e 2 /h at ip = tt/2. We find a 
Fano-suppressed Kondo plateau for G(e) similar to recent experiments. 
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Introduction. Recent measurements of transport 
through small semiconductor quantum dots have re- 
vealed interesting quantum-coherence phenomena such 
as the Kondo effect Jl]-^| in quantum dots strongly cou- 
pled to leads, Aharonov-Bohm (AB) oscillations || of the 
current through multiply-connected geometries contain- 
ing quantum dots, and Fano-type line shapes in multi- 
channel transport situations ||. In this letter, we study 
the combination of all these effects in an AB geometry 
which contains a spin-degenerate single-level quantum 
dot, see Fig. ||. 

Interference between resonant transport through the 
quantum dot and the direct channel gives rise to asym- 
metric line shapes in the linear conductance as a function 
of gate or bias voltage, the well-known Fano effect 0. 
Such line shapes have been observed recently in linear 
transport through multi-level quantum dots B, where 
the nature of the direct ( "reference" ) transmission path 
has not yet been fully clarified. Furthermore, scanning 
tunneling microscopy measurements of magnetic atoms 
on gold surfaces || yielded Fano lineshapes in the tun- 
neling density of states, which have been successfully ex- 
plained theoretically under the assumption that only 
conduction electrons participate in tunneling. The setup 
we propose (see Fig. [I]) has, however, the advantage that 
a controlled separate manipulation of both interfering 
paths is possible. Enhanced AB oscillations due to the 
Kondo effect at low temperature were predicted |l0| for 
a similar geometry. 

Here, we define two precise criteria for detecting Kondo 
correlations in a closed geometry. The first one is based 
on the interplay of Fano and Kondo physics. We find that 
the Fano line shape in the Kondo regime is qualitatively 
different as compared to a nonintcracting system, where 
the Kondo effect is absent. Furthermore, we find that 
the well-known Kondo plateau of increased conductance 
in the Coulomb-blockade regime is suppressed due to the 
Fano effect, and can even be inverted into a Kondo val- 
ley. We call this behaviour the Fano-Kondo effect. The 



second criterion addresses the scattering phase. In the 
Kondo regime, this phase is tt/2 associated with unitary 
transmission 0|. We predict that this can be detected 
in the AB phase <p = 2tt$ / $o dependence of the linear 
conductance G(ip), which shows a frequency doubling and 
a maximum at ip = tt/2. Here, $ is the enclosed flux and 
$o = hc/e denotes the flux quantum. This is remark- 
able since it demonstrates that the symmetry relation 
G(ip) = G(—ip) referred to as "phase locking", which is 
an exact property of (closed) two-terminal setups Egj , 
does not spoil the possibility to exhibit unitary scatter- 
ing in the Kondo-correlated state. Hence, it is not neces- 
sary to employ multi-terminal or open-geometry setups 
for that task, which have been studied in order to avoid 
phase locking. Their theoretical interpretation, however, 
relies on the assumption that multiple scattering events 
are absent, a condition which is difficult to control in ex- 
periments. This may be one reason for the discrepancy 
between recent phase shift measurements and theoret- 
ical predictions p3| . Another reason which favors the 
use of a closed-geometry setup from the practical point 
of view is the fact that the signal is strongly reduced in 
open geometries, since most of the emitted electrons are 
absorbed by terminals in the periphery. 
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FIG. 1. AB interferometer with a quantum dot embedded. 

The Model. The model we consider is depicted in 
Fig. [j]. Electrons driven through the device have either 
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to go through the quantum dot, described by an Ander- 
son model, or can be transfered directly. The area en- 
closed by the two paths is penetrated by a magnetic flux 
The Hamiltonian H = H L + Hr + H D + H T + H LR 
contains H r = Ylka e kr% crr (i'ka-r for the left and right 
lead, r = L/R. The isolated dot is described by 
Ho = e^jjiiff + C/n|?T,|, where n a = c\ r c l7l and e is 
the level energy measured from the Fermi energy of the 
leads. Tunneling between dot and leads is modeled by 
Ht = J2kcrr(^ r ral ljr c a + H.c), where we neglect the en- 
ergy dependence of the tunnel matrix elements V^/r. 
The intrinsic line width of the dot levels due to tun- 
nel coupling to the leads is (in the absence of the up- 
per arm) r = T L + T R with T L/R = 2ir\ V L/R \ 2 N L/R 
where N^/r is the density of states in the leads. The 
electron-electron interaction is accounted for by a charg- 
ing energy U — 2Eq for double occupancy. The trans- 
mission through the upper arm is modeled by the term 
H L r = Y,kqa( WelVa laR a q<?L+'H-£-)- We choose a gauge 
in which the AB flux $ enters via ip — 27r<&/ <&o only the 
tunnel matrix element for the direct transmission. 

General current formula. The current from the right 
lead is given by / = Ir = ed(fiR)/dt = i(e/h)([H ,tir}) . 
The latter expression yields Green's functions which in- 
volve Fermi operators for the leads and the dot, 



with T b = 4x/(l 



being the background trans- 
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with notations G< ka (t) = i{al aR a qaL (t)) and G% k>(r (t) = 

i(a^, cri? c (T (t)). The indices q and k label the states in the 
left and right lead, respectively. The index al indicates 
that a dot electron operator is involved (in our simple 
model there is only one dot level). The first (second) line 
of Eq. (0) describes electron transfer from the left (from 
the quantum dot) to the right lead or vice versa. This 
transfer can be a direct tunneling process or a complex 
trajectory through the entire device. 

Our goal is to derive a relation between the current 
and Green's functions involving dot operators only. To 
achieve this we employ the Keldysh technique for the 
Green's functions G < , G r , and G a , where G r and G a are 
the usual retarded and advanced Green's functions, re- 
spectively. We write down Dyson-like equations, collect 
all contributions (including those with multiple excur- 
sions of the electrons to the leads and the quantum dot 
and arbitrary high winding number around the enclosed 
flux), and make use of the conservation of total current. 
We get I = (2e/h) J duj T {uj)[f L {u) - f R (u)] with the 
total transmission probability T(u>) per spin given by 



T(lu) = T b + y/aT h R h cos (pTRe G r (u) 

~ [a(l-T b cosV) -T b ] flmG r M, (2) 



mission probability, R b = 1 - Tb, x = n 2 W 2 N L N R , 
and r = r/(l + x). Asymmetry in the coupling of 
the dot level to the left and right lead is parametrized 
by a = AT lT r/Y 2 (for resonant transmission the zero- 
temperature conductance through the quantum dot is 
G = ae 2 /h per spin). The first term in Eq. (Q) describes 
transmission through the upper arm. The second and 
third term represent both transport through the quan- 
tum dot as well as interference contributions. We em- 
phasize that the Green's function for the quantum dot 
entering Eq. (Q) has to be evaluated in the presence of 
the upper arm. The expansion of Eq. (|^) up to first order 
in r and W has been used (13] to address the effect of 
spin-flip processes on the suppression of the interference 
signal. We note that Eq. (||) can be generalized to even 
more complex scattering geometries within the formalism 
presented in Ref. fl5| . 

Our result can be viewed as a generalization of the 
Landauer-Biittiker formula to the interacting case. It 
includes interference effects and is especially suitable to 
describe linear-response transport. In this regime T(u>) 
is needed for zero bias voltage only, i.e., we can calcu- 
late the equilibrium Green's function G r by using the 
Numerical Renormalization Group (NRG) technique |l6| ] 
generalized to the case of two reservoir channels coupled 
to the dot O]. In the following, however, we focus on 
the T = case where, using Fermi liquid properties, all 
relevant information can be extracted from the dot occu- 
pation number calculated by NRG. 

Friedel sum rule. For {7 = 0, the Green's function 
can be calculated exactly. In equilibrium, the self-energy 
S(w), defined by G r (w) = l/[u}-e-E(w)], reads S(w) = 
— (T/2)^/a^/xcos (p — iT/2 and is independent of lj for 
a flat band in the reservoirs. In our calculations, this 
condition is satisfied due to D 3> T where D is the half 
bandwidth (for a discussion of narrow-band effects see 
Ref. |1^| ). For arbitrary U we use the Friedel sum rule 
M which at zero temperature yields ImE(O) = —T/2 
and a renormalized level position e + Re 11(0). The latter 
is related to the resonant scattering phase shift 6 rcs = 
7r(n)/2 of the dot (where (n) — J2a( n a) ^ s the total dot 
occupation) by 



e = p [e + Re E(0)] = cot S ICS . 



(3) 



With these properties we obtain the zero-temperature 
linear conductance G = 2ge 2 /h (or the dimensionless 
conductance g = T(0)) from Eq. (Q). By some algebra, 
it can be recast into the the generalized Fano form 



9 = Tx 



{e + qf 



■ 2 

sin ip 



with the Fano parameter 



q — cot S = \J aRb/Tb cos 



ip. 



(4) 



(5) 
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Here, we have introduced the nonresonant phase shift S 
which is due to scattering of free electrons at the weak 
link formed by the direct tunneling path. Eq. (|]) can 
then be expressed completely in terms of phase shifts: 

sin 2 (6 ros -5) . 2 . 2 

9 = T h r -2— 1- asm (psm d Tes . (6) 

sin o 

Discussion. Figure |^ shows the linear conductance as 
a function of level position e for different values of T, at 
vanishing AB phase <p = 0. (For all figures we choose 
symmetric coupling, a — I.) We find asymmetric line 
shapes of the peaks and dips around e = and e + U = 0, 
which indicates the presence of the Fano effect. 
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FIG. 2. Linear-response conductance at zero temperature 
as a function of level position e for different values of the 
background transmission Tb. We have used U/D = 0.5 and 
Y/D — 0.063, where the half bandwidth D — 1 has been 
taken as the unit of energy. Note that due to the small ratio 
Y/D, the results obtained here are cutoff-independent. The 
AB phase is if — 0. 

The Kondo effect shows up at Tb = as a large plateau 
of unitary transmission. At finite Tb the plateau survives 
but is reduced in height, similar to findings in recent ex- 
periments ||. However, in these experiments the nature 
of the direct tunneling path has not yet been clarified. In 
particular, the direct tunneling strength cannot be tuned 
independently from the level broadening, which makes a 
quantitative comparison difficult. It would, therefore, be 
highly desireable to realize the geometry shown in Fig. [l] 
in order to directly verify the reduction of the Kondo 
plateau. This reduction, which we call the Fano-Kondo 
effect, can be understood analytically from Eq. (Q. In 
the presence of interaction, the occupation of the dot is 
given by (n) w 1 in the whole Coulomb blockade regime 
— U + T < e < — T, i.e., the resonant phase shift is 
<5 rcs = tt/2, and Eq. (||) simplifies to 

35 roa=Tr /2 = a (l - T b cos 2 <p) , (7) 

or g = ai?b for ip = within the whole plateau. We 
emphasize that the Kondo effect does not break down. 



The reduction of the conductance plateau is rather due to 
interference between the lower (Kondo-) and the upper 
(reference) arm. For noninteracting systems, the phase 
is 7r/2 at resonance (e = 0) only, i.e., there is no plateau. 
The formation of a reduced Fano-Kondo plateau is there- 
fore a clear indication for Kondo correlations in the quan- 
tum dot under study. Far away from the plateau, we get 
fires — * 0, 7T, and thus g = Tb from Eq. ^ as expected. 
For arbitrary <5 res we get 

.9^=0 = [a+ (1 - a)T h ] sin 2 (<5 ros - S) . (8) 

This means that the conductance has a peak (of height 
2e 2 /h for a = 1 or Tb = 1) at fi Tes = 5 ± tt/2 and is zero 
for S Tes = S in agreement with Fig. ^. We emphasize the 
complementary behavior of the limiting cases of a closed 
(Tb = 0) and an open (Tb = 1) channel in the upper arm. 
For Tb = we obtain the usual Kondo plateau |l],|3| but 
for T> = 1 we have S = tt/2 and the conductance is zero 
in the Kondo valley but 2e 2 /h outside. 

We now turn to the dependence of the conductance 
on the AB phase ip and the question whether the phase 
8 les = 7r/2 for unitary transmission in the Kondo regime 
can be detected from the AB oscillations of the linear 
conductance for a closed AB interferometer. We have 
calculated the shape of the AB oscillations for different 
gate voltages in and outside the Kondo regime, as shown 
in Fig. ||. Because of the phase locking property G(ip) — 
G(—(p), we have plotted only ip 6 [0, n]. The inset shows 
the lowest coefficients of the Fourier expansion G(<p) — 
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FIG. 3. AB oscillations of the linear conductance for dif- 
ferent level positions e. Parameters are identical to those in 
Fig. |i| with a background transmission of Tb = 0.3. Inset: ab- 
solute value of Fourier coefficients, \a n \. In the Kondo regime, 
the first harmonic vanishes and a-z dominates, thus effectively 
doubling the frequency of the AB oscillations. 

Outside the Kondo regime (e/D = and —0.5), the 
AB oscillations are dominated by the lowest harmonic 
(period 2tt), with global extrema at ip = 0, tt. A drastic 
change occurs as soon as the dot is tuned into the Kondo 
regime: then, the AB oscillations show a maximum at 
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if > = 7r/2 with the universal conductance 2ae 2 /h inde- 
pendent of the value of the background transmission Tb. 
As a consequence (see inset of Fig. ||), the n = 1 Fourier 
coefficient vanishes, and the second harmonic dominates, 
thus effectively doubling the oscillation frequency. Due to 
Kondo screening, this frequency doubling persists over a 
finite range of gate voltages corresponding to the con- 
ductance plateau. This effect, therefore, yields a conve- 
niently measurable and precise criterion for Kondo cor- 
relations in a quantum dot. 

Figure ^illustrates the "pinning" of the AB maximum. 
We show G(e) for different AB phases ip. For ip = tt/2 
we recover a Kondo plateau with height 



9ip=Tr/2 



a + (T b - a) cos 2 5 T , 



(9) 



A special case is Tb = a where the conductance is given 
by 2ae 2 /h for all level positions. Furthermore we note 
the property G(e, ip) = G{— e — U, f + n) which follows 
from particle-hole symmetry. 
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FIG. 4. Dependence of the Fano lineshape on the AB phase 
for the same parameters as in Fig. ^ and a background trans- 
mission of Tb = 0.3. 

Summary. We have studied the interplay of Kondo 
and Fano physics in the most basic model describing a 
closed AB interferometer which contains an interacting 
quantum dot in one of its arms. We derived a current 
formula and used Friedel sum rule to calculate the ex- 
act linear conductance at zero temperature with the help 
of the Numerical Renormalization Group generalized to 
two-channel systems. We found a characteristic Fano- 
suppressed Kondo plateau. Furthermore, we demon- 
strated that unitary transmission of a Kondo-correlated 
state can be identified by a frequency doubling of the 
AB oscillations and a conductance maximum at ip = n/2 
with universal height. 
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Note. While this paper was written up, a preprint |ll]] 
was published, in which a similar formula as Eq. (|^) has 
been presented for the same model. However, we believe 
that the result of Ref. jl9| is incorrect since it differs from 
our result by a factor (l+x) in the second and third term. 
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